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^Jj Abstract 

A quantum cosmological model with radiation and a dilaton scalar field 
is analysed. The Wheeler-deWitt equation in the mini-superspace induces a 
Schrodinger equation, which can be solved. An explicit wavepacket is con- 
structed for a particular choice of the ordering factor. A consistent solution 
is possible only when the scalar field is a phantom field. Moreover, although 
the wavepacket is time dependent, a Bohmian analysis allows to extract a 
bouncing behaviour for the scale factor. 

1 Introduction 

A scalar field is the simplest fundamental matter field that we may introduce 
in a cosmological model, since it adds essentially just one degree of freedom 
and it is naturally invariant under coordinate transformation. Even though, 
we can quote two different couplings of a scalar field: the minimal coupling, 
leading to the so-called Einstein's frame; the non-minimal coupling, corre- 
sponding to the Jordan's frame. Of course, the minimal coupling is the most 
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simple one. But, the non-minimal coupling is one of the predictions of string 
theory: the scalar field emerging from this theory, called dilaton field, appears 
non- minimally coupled to the gravity term pQ. Moreover, being an interest- 
ing alternative to the General Relativity theory, the Brans-Dicke theory [2], 
is based on the non-minimal coupling. 

Both frames, the Einstein and Jordan ones, may be connected by a confor- 
mal transformation. There is an intensive discussion in the literature about 
the meaning of such connection, in the sense that it may be just a mathe- 
matical mapping of one system into another, or may hide a deep physical 
meaning, see for example [3] and references therein. For some observables, it 
is clear that the predictions in one frame are not equivalent to the predictions 
in the other frame. In quantum cosmology this problem has been treated in 
references jH [5] , showing that it is possible - even at quantum level - to mapp 
the equations in one frame into another frame, but with different predictions 
to the evolution of the universe. 

In references [U |5] a system consisted of gravity and a scalar field, mini- 
mally or non-minimally coupled has been considered. In this case, the anal- 
ysis present a major challenge, typical of quantum cosmology: how to obtain 
predictions for the evolution of the universe as function of time? In fact, it is 
well known that, considering the Einstein-Hilbert action, in quantizing it the 
resulting equation, the Wheeler-de Witt equation, has no time parameter due 
to the invariance of the gravitational system under time reparametrization. 
The introduction of a scalar field does not change this situation. In those 
references, the WKB method has been used in order to obtain predictions 
for the evolution of the universe from the wavefunction determined from the 
(timeless) Wheeler-de Witt in the minisuperspace. Concerning this question, 
see reference jS] and references therein. 

Another possibility to obtain the time evolution of a quantum cosmo- 
logical system in the minisuperspace is to introduce matter in the form of 
a fluid, employing the Schutz description based in terms of potential that 
conveys the degrees of freedom of the fluid [9j [10] . It has been shown that 
such description leads, always in the minisuperspace, to a Schrodinger-like 
equation, with the degrees of freedom of the matter playing the role of time. 
This proposal has been presented in reference [11], and extensively analysed 
in reference [T2] . 

Our aim in this work is to address the problem of time evolution of the 
universe in the presence of a dilaton-like field from a quantum cosmologi- 
cal perspective. Since we are interested mainly in the primordial universe, 
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a radiative fluid will be introduced. This may allow to recover the time 
variable through the employment of the Schutz formulation. Performing a 
conformal transformation (which does not affect the radiative fluid, since it 
is conformal invariant), we re- write the dilaton-gravity system in the Ein- 
stein's frame, from which the Wheeler-de Witt equation in the minisuper- 
space is constructed, resulting in a Schrodinger-like equation. From this 
Schrodinger-like equation, an explicit solution is obtained by using a specific 
ordering factor. A wavepacket is determined, but its norm is time- dependent. 
Hence, it does not fit in the usual many-worlds interpretation of quantum 
mechanics [7J |6], but it admits a Bohmian analysis [6j Ej. The Bohmian 
trajectories contain universes that are singurality-free. 

This paper is organized as follows. In the next section, we quantize 
the dilaton-gravity-radiation system in the mini-superspace. In section 3, a 
wavepacket is constructed and a (formal) many-world analysis is performed. 
In section 4 the Bohmian trajectories are studied. In section 5 we present 
our conclusions. 

2 Dilaton-gravity system with radiative fluid 

Let us consider the non-minimal coupling of gravity and a scalar field, rep- 
resented by the Brans-Dicke theory: 



where L m is the matter Lagragian supposed to be conformal invariant (a ra- 
diative fluid). In the strict string dilatonic case we have u = —1. This defines 
the theory in the Jordan's frame. Performing a conformal transformation 
such that = we transpose the action ([TJ to the corresponding 

expression written in the Einstein's frame: 



where to = cu + 3/2. 

Restricting ourselves to the FLRW metric, in an appropriate coordinate 
system, the line element can be written as 




(1) 





ds 2 = N(t) 2 dt 2 - a(t) 2 'j ij dx i dx j 



(3) 
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where N(t) is the lapse function, a(t) is the scale factor and 7^ is the in- 
duced metric of the homogeneous and isotropic spatial hypersurfaces with 
curvature k — 0, ±1. From now on, we will fix k = 0. With this metric, the 
gravitational Lagrangian becomes, 



r V a- 



\ 2 

a 1 a\ a N 



a \a J a N 



l2 



"J*} > ( 4 ) 



where Vq is a constant and can be interpreted as the physical volume of 
the compact universe divided by a 3 . Since we shall have an identical multi- 
plicative constant in front of the matter lagrangian we can drop it from our 
analysis (this can also be understood as a normalization of the fields). Hence, 
discarding a surface term, the gravitational Lagrangian can be written as, 
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C G = -| 6 aa 2 - uja^j. (5) 
a = y / jwjln0, (6) 



k2 



Defining, 
we obtain, 



where e = ±1 according u is positive (upper sign) or negative (lower sign). 
The canonical momenta associated with the scale factor and the scalar field 
are respectively: 

Pa = 12- , V, = ~2e— . (8) 
This leads to the following expression in terms of the conjugate momentum: 

£™+»*-tf{l£-4} . (9) 

Considering a radiative matter component (for the computation of the con- 
jugate momentum associated with the fluid, see references [TTIIT2] ). the total 
Hamiltonian is: 

2f = Jv/i^-c4--V (10) 



24 a 4a 3 
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The resulting Schrodinger equation is, 

where we made the redefinition — > a and ^ — > T. 

Two questions related to ( ITTj) appear immediately. First, what is the sign 
of e? If it is positive, we have a hyperbolic Schrodinger equation. This implies 
that the "energy" E is not positive defined, as it can be seen by multiplying 
the Schrodinger equation by and integrating in a and a. Moreover, in this 
case, the argument of the Bessel function can become imaginary, what may 
not be very serious, but may pose some problems to construct the wavepacket, 
since the integration in E must be done along all real axis, and for E < the 
Bessel function becomes modified Bessel function. On the other hand, if the 
sign is —1 the positivity of E seems to be assured, and the construction of the 
wavepacket seems not problematic, since the integration on E is done only 
in the positive semi-axis. For all these reasons, for the moment, we consider 
the case e = — 1. It corresponds to a phantom scalar behaviour induced to 
the dilaton field. 

Second, we have also not taken into account a possible factor ordering. 
If such factor ordering is introduced, with e = —1, we have the following 
Schrodinger equation: 

d 2 m pd^ 1 d 2 ^ .dm 

Even though there is no unique way to chose the factor ordering, a possible 
consistent and adequate choice is the covariant ordering which is invariant 
through fields re-definitions. The minisuperspace for our model is two di- 
mensional. Thus, using the covariant factor ordering, we can re-write the 
Wheeler-DeWitt equation in a Klein-Gordon like equation by defining the 
minisuperspace metric such that eq. (fl2l) reads 

^MvW^,)* = -^ , (13) 

where /i, v = 0, 1 and the coordinates represents the fields, i.e. x° = a and 
x 1 = a. The only possible value for p in eq. (I12j) which allows the covariant 
factor ordering is p — 1. In this case, we have 

9lM/= \0 a 2 ) ^ V9 = a (14) 
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and the Wheeler-DeWitt equation reads 

"1 d ( d \ 1 d 2 1 T M 
_ad-a[ a d-a) + ^\ 'df ' (15) 

This equation can be solved using the technique of separation of variables. 
Hence, through the ansatz 

^(a,a,T) = <j)(a)e ika e- iET , (16) 

we obtain, 

0" + -+<U-^)0 = O, (17) 

where the prime means derivative with respect to a. Noticing that this is 
just a Bessel's equation, its solution reads 

* = AJ u (VEa)e ika e- iET , v = k , (18) 

with A being a normalization constant. 



3 Wavepacket and expectation values 

Let us construct a wavepacket by choosing conveniently the function A 
A(E, k). We have, for (PJ), 



OO POO 



*(a,cr,T)= I / r u+1 e-^ +lT)r e- {a - iu)k J k {ra)dkdr, (19) 
'o Jo 



where r = yE and a, 7 are postive constants. Using the formula (6.631-4) 
from [T3] , the final result is: 



®(a,a,T) = C 



e *Bm 



B(T) g a (a,B,a) 



(20) 



where 



B(T) = (j + iT), g a (a, a,T) = -a + In 



2B(T) 



±ia, (21) 
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and C is a normalization constant. Remark that, in order to give physical 
meaning for this wave packet, the condition $lg a (a,a,T) < must be satis- 
fied, assuring that the integral in the separation parameter k is convergent. 
Let us calculate this normalization constant. Unitarity requires, 



N — ^*^>dada = 1. (22) 



But, 



Remark that, 



where 



= c 2 . (23) 

2B*B g*g v ; 



B = 7 + %T = De ie , (24) 



D = vV + T 2 = VB^B, (25) 
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6 = arctan(-). (26) 



Hence, 



g* g = h 2 + (a-6)\ h = -a + \n[- - T ). (27) 

For the scalar field a and for a we must compute a double integral. In doing 
so, we must fix the range of a. In principle, it is from — oo < a < +oo, 
implying < <p < +oo, that is, a positive gravitational coupling. The norm 
of the wave function becomes 

M _ C 2 f°° f +co e^^F da 

N ~ (2BB*y/ 2 J (2B*B)l[h?+{a-9?\ da - (28) 



Under the substitutions, 



m = r, v = -9 + a, (29) 

(2B*5)^ 



the integral above becomes 



Q2 poo r+oo — 7U 2 

N = ru d *m/2 / / TTT-^ dudv - ( 3 °) 
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But 



+00 



dv 71 



-00 



h 2 + v 2 h' 



(31) 



N = WW^" Jo ; ( u \ dU = TBBW^ ^ 



Hence, 

C 2 f°° e-T* C 2 



where I\ is the definite integral, 



a + ]n( f 



e 



Ji = / ^rfn. (33) 

a + lnff 

The wavepacket determined above is time dependent and does not convey 
unitarity. 

As a formal exercise, even in absence of unitarity we evaluate the expec- 
tation value of the scalar field and of the scale factor for this wavepacket: 

roc r+°°y* a y dada 

<a>r = Vr- ' (34) 

JO J — 00 

C II™ H" a Vila da 
<a >T = J-r^dada ■ <35) 

JO J ~ 00 

Let us first evaluate the expectation value for the scalar field. In the numer- 
ator of (1531). we have: 



a da f +co (v + 9)dv 



h 2 + ( a _g)2 J ^ h 2 + v 2 

The first integral is zero, and the second one leads to, 



(36) 



+ °° _ /•+" dv _ rr 

-oo h 2 + ia _ d ) 2 h 2 + v 2 h' 1 1 

Re-inserting in the expression for the scale factor, we find: 

< a >t= arctan( — ] . (38) 



\7 
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For the scale factor, integrating in a we have: 

f°° ae- 2 ^^ 

<a>T=7, l T,^B- da - (39) 

Performing, as before, the substitution 

n n 

u, (40) 



we obtain: 



where 



Hence, 



< a > T = V7 2 + T 2 \e~ 2 ^du = h^/^ + T^ 
Jo " h 



u 



(41) 



^e-^du. (42) 



<a> T = a ( 7 2 + T 2 ) 1/2 , (43) 

which is formally the same solution as when there is no scalar field [T2] . 
Remark that the norm of the wavepacket is time dependent and a unitary 
framework can not be invoked. 

The fact that the result for the scale factor is the same, essentially, as 
when the scalar field is absent may be seen as contradictory. But, we remark 
that in the Friedmann equation it appears the logarithmic derivative of the 
scale factor and the derivative of the scalar field. Hence, the presence of the 
scale field just change the value of a with respect to the case it is absent. 



4 Bohmian trajectories 

In order to extract some physical insight on the model developed until here, 
given the nature of the wavepacket employed in the previous section, we 
will use a Bohmian analysis [6]. Note that the computation of the Bohmian 
trajectories makes sense even in absence of unitarity, this being one of the 
main features of the ontological interpretation of quantum mechanics. 

The Bohm-de Broglie interpretation is most easily understood through 
the polar form of the wave function. Indeed, decomposing the wave function 
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as \& = Re lS , the Wheeler-DeWitt equation (|T5"|) split in two real non-linear 
equations for the two real function R(a, a, T) and S(a, a, T) 



dR 2 
dT 

dS_ 
df 



1 9 



dS 



-— ( B? 2a— 

oa 



a da 

fosV 

\da) 



1 d 



a da 



2dS 



-— Rf--^- 



a da 



ldS\- 
a da J 



with 



i__9_ 

a 9a 



9/2 

9a 



+ 



1 9 2 fl 
a 2 9a 2 



(44) 
(45) 

(46) 



Equation ( 145]) is a modified Hamilton- Jacobi equation with the presence 
of the quantum potential Q. Through this equation, we can identify the 
momenta as 



P -™ 
da 



P -™ 
a ~ da 



(47) 



The metric in the minisuperspace (fl4|) allows to define covariant derivatives, 
see reference [13]. Therefore, equation ( 144|) represents a continuity equation 
written in that minisuperspace, with R 2 playing the role of a probability den- 
sity and we can directly read the associated velocity within each parenthesis, 
i.e. 
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[XV 



dS 



as 



dx u da 
The wavefunction can be written as, 



a 



2_as 

a 2 da 



V(a,a,T) = A 



e 4s 



a z B* 

e iBB* 
A ^^g* B* 



gB g*gB*B 
f(a,a,T)e^ T \ 



(48) 



with, 



/(a,<7, T) 
S(a, a, T) 



A- 



g iBB* 

\Jg*gB*B' 
a 2 T 
4( 7 2 + T 2 ) 



arctan — + arctan 
\7/ 



( 



-d + a* 



\ h 



(49) 
(50) 
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Remembering that, 



h = —a + \n[ = 9ft qr(a, &)■ (51) 

and that, for a radiative fluid N = a in (JHJ) [12], taking into account the time 
and field redefinitions, we have the following expressions for the Bohmian 
trajectories: 

aT 2 9 -a 

+ ~ ,2 , ta ,v, . ( 5 2) 



( 7 2 + T 2 ) a /i 2 + (0 - ^ 2 
2 

" /i 2 + (0 - tr 



* = —7^^ 35 » ( 53 ) 



where and /i are given by (I26|l5ip . These equations must be integrated 
numerically. 

The numerical solutions of f)52|53p provide an interesting perspective de- 
pending on the initial conditions and on the free parameters of the model, 
namely 7 and a. We can show that there is at least a class of non-singular so- 
lutions, satisfying the wavepacket conditions of consistency. More precisely, 
in figure 1, the a non-singular solution is displayed. The function h is also 
plotted, showing that the condition for the convergence of the wavepacket is 
satisfied. In the same figure we plot the results obtained use the expectation 
value of the previous section, and the numerical solutions for the Bohmian 
trajectories. The Bohmian trajectories predict an asymmetric bounce, while 
the expectation value display a perfect symmetric bounce. 



5 Conclusions 

The possibility to make predictions and subsequently extract a family of solu- 
tions (trajectories in the minisuperspace) is a recurrent problem in quantum 
cosmology. Although the well-known WKB method, including decoherence 
features, has provided a vast range of particular results, there is an alternative 

6 It can already be stated that the solutions are not equivalent to the expectation values 
found before, a consequence of the absence of unitarity: when the wavepacket is unitary, 
the Bohmian trajectories reproduces the expectation values. 
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Figure 1: In the left panel, the behaviour of the scale factor for 7 = 100 and 
a = 50 is shown. In the central panel, the function h is displayed. In the 
right panel, the solution obtained using the Bohmian trajectories (blue) is 
compared with the solution using the expectation value (red). 

framework by means of transforming the WdW equation into a Schrodinger 
equation, with the time variable induced by a matter component [HI EU ED]- 
Quite interesting results have appeared in the literature [12] using this pro- 
cedure. 

In this paper, we investigated how we could obtain some predictions with 
that alternative framework, using a system characterized by a dilaton field 
and radiation expressed using the Schutz variables. The dynamics of the 
radiative fluid implies a time variable. 

We found that, in order the Schrodinger equation to be elliptic, leading to 
a positive energy spectrum, the dilaton field must have a phantom behaviour. 
On the other hand, the construction of a quasi gaussian superposition, con- 
jugated with some convergence criteria, lead to a wavepacket with a time 
dependent norm and hence unitarity could not be invoked. Nevertheless, 
this still allows us to investigate this wavepacket under a Bohmian perspec- 
tive. Interestingly, we found a bouncing behaviour for the scale factor, i.e., 
the singularity is avoided. In this sense we extend the results of reference 
[15], where the imposition of unitarity led to a conclusion that the dilaton 
field should be constant. At same time, such approach open new perspective 
concerning anisotropic quantum models |16j . 

There are still many open directions to explore in this program. In partic- 
ular, we must look for, (i) other wavepackets, presumably by using numerical 
methods, (ii) study the dynamics in the Jordan frame. We hope to address 
these problems in a future work. 
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